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Abstract 

When the crew exploration vehicle (CEV) is launched, the spacecraft adaptor (SA) fairings that cover 
the CEV service module (SM) are exposed to aero heating. Thermal analysis is performed to compute the 
fairing temperatures and to investigate whether the temperatures are within the material limits for nominal 
ascent aero heating case. Heating rates from Thermal Environment 3 (TE3) aero heating analysis 
computed by engineers at Marshall Space Flight Center (MSFC) are used in the thermal analysis. Both 
MSC Patran 2007rlb/Pthermal and C&R Thermal Desktop 5.1/Sinda models are built to validate each 
other. The numerical results are also compared with those reported by Lockheed Martin (LM) and show a 
reasonably good agreement. 


Introduction 

The ascent heating on the crew exploration vehicle (CEV) is analyzed by using computational fluid 
dynamics (CFD) and engineering codes at MSFC. Thermal Environment 3 (TE3) heating data is used in 
this work. One of the major concerns is with the SA fairings covering the CEV service module (SM) and 
the SM/crew launch vehicle (CLV) flange interface. The TE3 heating rate is a function of time, wall 
temperature, and the spatial locations. Two commercial software packages, Thermal Desktop (TD) 5.1 
and MSC Patran 2007rlb, are used in this analysis. TD/Sinda has been widely used for spacecraft-related 
thermal analysis including orbital heating. TD is the pre- and post-processor for Sinda, which is a finite- 
difference-based solver. In TD, the geometry is built and meshed, the boundary conditions are defined, 
and then Sinda is used to compute temperatures. To define the BCs, a conductor can be defined using the 
node-to-node, node-to-nodes, or node-to-surface options. A contactor can be defined using the surface-to- 
surface or edge-to-surface options. Heat loads can be defined on nodes, surfaces or solids. MSC Pthermal 
is a finite element-based thermal solver. MSC Patran is the pre- and post-processor for Pthermal. It can 
import geometry from ProE parts or assembly files, IGES files, and step files. Regarding the BCs, the 
convection, contact resistance, and heat load can be imposed in different ways. Orbital heating can not be 
modeled in MSC Pthermal. 

The TE3 heating rate will be presented first and the simplification of the data representation is 
described. Then the thermal models and results obtained by using TD/Sinda, and MSC Patran/pthermal 
are presented. The temperature results are compared; finally, the conclusion is given. 

The TE3 Ascent Heating Rate on SA 

The ALAS-1 1 TE3 ascent heating rates are gathered at certain body points (BPs) on CEV and CLV. 
On the section of the SA, there are three BPs along the axial direction, denoted as BP 016x03, 016x05, 
016x09, and are shown in figure 1(a). On the SM/CLV flange interface, the BPs are denoted as 021x20, 
021x21, 021x22, 021x23, and 021x24 and are shown in figure 1(b), where x = 0,1, 2, ...7 for eight points 
uniformly distributed along circumferential direction. In figure 1(c), the geometry of the SM/CLV 
interface is shown and each part name will be used in the following discussion. 
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Figure 1 . — The TE3 BPs on SA and SM/CLV interface flange, and the name of each part on SM/CLV interface 
for future reference, (a) BPs on SA. (b) SM/CLV interface BPs. (c) SM/CLV interface part names. 


As an example, the heating rates at all BPs on the SA are plotted in figure 2 for the cold wall case. 

The heating rate on the SA is nearly uniform along the axial and circumferential directions. To simplify 
the thermal analysis, the heating rate at BP 016703 (x = 802 in., 0 = 315°) has the highest heating rate for 
most of the time, and will be used over the entire surface of the outer face sheet (OFS) of the SA fairings. 
For the BPs on the SM/CLV interface, the detailed heating rate is referred to (ref. 1) and will not be 
plotted here. The heating rate at 0 = 0° is used for all different axial locations. The heating rates provided 
in TE3 data include gas convection, gas radiation, particle convection, and particle kinetic energy heat 
transfer. For the ascent heating on the SA, only the gas convection contributes to the heating rate. 
Therefore, the total heating rate can be defined as 

9 total — (^rec — ^wall ) ( 1 ) 


where h c is the enthalpy -based heat transfer coefficient, H rcc is the gas recovery enthalpy, // wa ii is the 
gaseous wall enthalpy and defined as 

7f wall = 0.2345r wall +9.786xl0- 6 r4 n + 943.6/r wall -1.57 (2) 
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Figure 2. — TE3 heating rate at the body points on SA for the cold wall case, (a) BP 016x03. (b) BP 016x05. 

(c) BP 016x09. 

with r wa n being the wall temperature in Rankine. The given time history of h c changes slightly when the 
wall temperature changes from 0 to 2000 °F. It is obvious that g total is not a linear function of 7" wa ii- With 
the given ascent heating data, g total is plotted as the function of 7 wa n in figure 3. It can be seen that g total 
is almost a linear function of r wa n when r wa n<760 °F. Since the material temperature limit on the SA is far 
below 760 °F, it can be assumed that 


^wall - 0.2345r wall (3) 

Based on this assumption, the enthalpy-based heat transfer coefficient (h c ) can be easily converted into 
the temperature-based heat transfer coefficient (h t ) as follows: 

^ t ~ c ph c -> ^rec — ^rec / c p •> ^wall — ^wall / c p (4) 

where c p = 0.2345 Btu/lbm-°F for air, and T XQC is the air recovery temperature. The so-obtained h t and T rec 
are plotted in figure 4. It shows that h t changes when the wall temperature is different. The heat rates 
computed based on /z,and/z c are plotted in figure 5, which shows that h t gives slightly higher heat flux 
than h c . The convection BC on the fairing surface is only time dependent when using h t . Thermal analysis 
presented here is on the conservative side. It is the worst case for the temperatures on the SA fairings. 


NAS A/TM— 2009-2 1 5474 


3 


h, Btu/ft 2 -s-F, temperature based 


0.7 


cp 

C\! 

£ 


0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

0.0 


_L 


_L 







j j 








1 I 








L i 








7 




















L I 

_"j 







■ i 

i 



_L 


_L 


_L 


_L 


_L 


_L 


_L 


200 400 600 


800 1000 1200 1400 1600 1800 2000 

7"wall- 


Figure 3. — The heating rate versus the wall temperature at BP 016709 at t = 100 s. 
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Figure 5. — The comparison of the enthalpy-based and temperature- 
based heating rate. 
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Thermal Model in Thermal Desktop 


In the TD model, the ACIS file is imported to build TD surfaces and solids. The SA fairings are 
formed by an OFS, a Honeycomb (H/C) core, and an inner face sheet (IFS). Each has three pieces along 
the circumferential direction. The OFS and IFS are 0.0424-in.- thick composite, and H/C core is 1.5-in.- 
thick aluminum. Each piece is modeled as a cylindrical surface with a defined thickness. 

The convection BCs are defined using the node-to-surface option. The air recovery temperature is 
defined at the node, and h t is defined on the OFS of SA fairing and SM/CLV interface. The contact 
resistance between the inner surface of the OFS to H/C core, and the outer surface of the IFS to H/C core 
is h = 0.694 Btu/hr-in 2 -°F for low resistance. The IFS to SM ring (edge to surface) is h = 0.0014 Btu/hr- 
in-°F, assuming that a minimum heat transfers to the ring from IFS and the overlap between IFS and ring 
is 2.0 in. The SM ring to back cone (edge to surface) is h = 14 Btu/hr-in-°F for low resistance. The flange 
to ring on CLV, flange to ring on SM, and CLV flange to SM flange (surface to surface) is h = 6.94 
Btu/hr-in 2 -°F for low resistance. The outer surface of the SA fairing radiates to the ambient air at 
T=50 °F, and the inner surface of the fairing radiates to the radiator at T= 70 °F. 

The SA fairings jettison at t = 150 s. The temperature contour at t = 100.8 s is plotted in figure 6 for 
the OFS, H/C core, and IFS, when the fairings reach the maximum temperature. Figure 7 shows the 
corresponding contour for the SM/CLV flange interface at t = 172.8 s. The time history of the temperature 
on OFS, H/C core, IFS, and SM/CLV interface is plotted in figure 8, showing that the temperature could 
reach 275 °F on the OFS and 225 °F on the SM ring. Note that no heat load is defined on the hinge 
(fig. 7). Since the hinge is a protuberance above the SM ring, the local area next to the hinge could have a 
higher heating rate. 
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Figure 6— TD results of OFS, H/C core, and IFS. Figure 7— TD result of SM ring and SM/CLV flange interface. 



Figure 8. — Time history of the temperature on OFS, H/C, IFS, SM ring, and SM/CLV flange interface (TD results). 
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Thermal Model in MSC Patran/Pthermal 


In the MSC Patran model, there are 1 16,552 nodes and 290,249 elements. The BCs are the same as 
those used in the TD model. The Pthermal temperature contour on the SA fairings, IFS, H/C core, and 
OFS at t = 100 s are plotted in figure 9. The MSC Patran model has more accurate modeling of the actual 
geometry of OFS, H/C core, and IFS. It shows more temperature variations across each part since it has 
different thicknesses even with the same ascent heating rates are applied. The temperature time history on 
the OFS, H/C core, and IFS are plotted in figure 10. It shows that the maximum temperature on the 
fairings reaches 300 °F at t = 100 s. The corresponding results for the SM/CLV interface structures are 
plotted in figures 1 1 and 12, showing that the maximum temperature there reaches 230 °F at t = 170 s. 
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Figure 9. — MSC Patran result of OFS, H/C, and IFS at t = 100 s. 
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Figure 10. — Time history of the temperature on OFS, IFS, and H/C (MSC Patran results). 
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Figure 11. — MSC Patran/pthermal result of SM ring and SM/CLV flange interface. 


The detailed temperature comparison between the TD results, MSC Patran results, and LM reported 
results (ref. 2) are shown in table I, showing a reasonable agreement. The maximum allowable 
temperature is also listed. For the SA fairings, the temperature is within the limit. However, the 
temperatures on SM/CLV flanges and SM ring are close to or above the temperature limit. Although the 
analysis presented here used a conservative heating rate, the hinge effect is not included. The hinge is a 
large protuberance on the SA fairing and will results higher heat transfer. Necessary thermal protections 
will be considered at the local area on the SM/CLV flanges and SM ring close to the hinge. Further 
analysis on the thermal protection is necessary. 
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Figure 12. — MSC Patran/pthermal result of SM ring and SM/CLV flange interface. 


TABLE I.— COMPARISON OF MAXIMUM TEMPERATURE ON SA FAIRING 
AND SM/CLV FLANGE INTERFACE 



Max temperature a , 
°F 

Max allowed 
temperature, 
°F 


LM result 

Current 
TD result 

Current MSC 
Patran result 


Fairing (OFS) 

272 

280 

280 

325 to 400 

SM/CLV flange 

184 

186 

196 

150 to 180 

SM ring 

200 

225 

232 

220 

LSC 

200 

NA 

210 

220 


a No hinge localized heating (with margin of 1.35 included in TE3 data) 


Conclusions 

The ascent heating thermal analysis on SA fairings and SM/CLV flange interface are performed using 
MSC Patran and C&R TD. Both TD and MSC Patran results agree reasonably well with those reported by 
LM on fairings and SM/CLV flange interface. Both results predict higher temperature than LM reported 
because of the use of temperature -based heat transfer coefficient in the analysis, which gives higher heat 
flux. The SA fairing reaches maximum temperature (280 °F) at / = 100 s. The hottest spots are next to 
hinges, and the areas that are not covered by honeycomb (as shown in MSC Patran model, 320 °F). The 
SM flange and ring and LSC reach maximum temperature at / = 170 s. The temperature at the comer of 
the ring could reach 253 °F. With the hinge effects, the fairing temperature under the hinges could get 
higher than predicted here. Necessary thermal protection is needed for the areas on the SM/CLV flanges 
and SM ring that are next to the hinge. 
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